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ABSTRACT 



Context. Absorption of ultraviolet radiation by water ice coating interstellar grains can lead to dissociation and desorption of the ice 
molecules. These processes are thought to be important in the gas-grain chemistry in molecular clouds and protoplanetary disks, but 
very few quantitative studies exist. 

Aims. We compute the photodesorption efficiencies of amorphous water ice and elucidate the mechanisms by which desorption oc- 
curs. 

Methods. Classical molecular dynamics calculations were performed for a compact amorphous ice surface at 10 K thought to be 
representative of interstellar ice. Dissociation and desorption of H2O molecules in the top six monolayers are considered following 
absorption into the first excited electronic state with photons in the 1300-1500 A range. The trajectories of the H and OH photofrag- 
ments are followed until they escape or become trapped in the ice. 

Results. The probability for H2O desorption per absorbed UV photon is 0.5-1% in the top three monolayers, then decreases to 0.03% 
in the next two monolayers, and is negligible deeper into the ice. The main H2O removal mechanism in the top two monolayers is 
through separate desorption of H and OH fragments. Removal of H2O molecules from the ice, either as H2O itself or its products, 
has a total probability of 2-3% per absorbed UV photon in the top two monolayers. In the third monolayer the probability is about 
1% and deeper into the ice the probability of photodesorption falling to insignificant numbers. The probability of any removal of H2O 
per incident photon is estimated to be 3.7 x 10""^, with the probability for photodesorption of intact H2O molecules being 1.4 x 10~^ 
per incident photon. When no desorption occurs, the H and OH products can travel up to 70 and 60 A inside or on top of the surface, 
respectively, during which they can react with other species, such as CO, before they become trapped. 

Key words, astrochemistry - molecular data - ISM: molecules 



1. Introduction 

Ices are a major reservoir of the heavy elements in a variety 
of astrophysical envir onments, rang ing from cold and dense 
molecular c l ouds ( e .g.. IWillner et alJ \V 
Smith et alJ 1 19891: |Murakawa et al " ^ 
2004 and proto planetary disks 



1982; Whi ttet et al 
^000; PontoDPidan et al. 
(Pontoppidan et al. 2005; 



Terada et al o the icy bodies in our own solar system such 

as comets (e .g. jMu mma et aD ll993') and Kuiper Belt Objects 
(jjewitt & Lu u 2004). In star-forming clouds, the fraction of car- 
bon and oxygen locked up in ice is comparable to that in gaseous 
CO (IPontoppidan 2006; van Dishoeck et al. 1996), whereas at 
the centers of cold pre- stellar co res more than 90% of the heavy 
eleme nts can be frozen out (e.g. JCaselli et al .111 9991: iBer gin et al.l 
I2OO2I) . Similarly, the cold midplanes of protoplanetary disks 
around young stars are largely devoid of gaseous molecules 
other than H2, H3 and their isoto pologues (e.g., Aikawa et al. 



l2002l:ICeccarelli & Dominikll2005l) . Thus, a good understanding 
of how molecules adsorb and desorb from the grains is critical to 
describe the chemistry in regions in which stars and planets are 
forming. 

The importance of ultraviolet (UV) radiation in aflTecting 
interstellar ices is heavily debated in the literature. On the 
one hand, the large extinctions of 100 mag or more along the 
lines of sight where ices are detected prevent UV radiation 



from penetrating deep into the clouds (e.g., lEhrenfreund et al.l 
[2001: Stauber et al. 2004), unless there are cavities throug h 
which the stellar UV photons can escape (ISpaans etal.lll995l) . 
Thus, the bulk of the ices are thought to be shielded from 
both external and internal radiation sources in whic h case pho- 
todesorption is thought to be unimportant (e.g., 'Leger et aP 
1985; Hartquist & Williams 1990). On the other hand, UV pho- 
tons are produced locally throughout the cloud by the interac- 
tion of cosmic rays with the gas, albeit at a level about 10^ 
times less than that of the general inters tellar radiation field 
(iPrasad & Tarafdaj l 19831: Ishen et al.l l2004). In addition. X-rays 
from young stars penetrate much further into their surroundings 
than UV and can produce loc a l UV photons throu gh a similar 
process ("Dalgarn o et al.l 119991: IStauber et al.l 12005 ). Moreover, 
the observed emission of optically thick millimeter lines from 
gaseous molecules is often dominated by the outer layers of the 
cloud where UV photons play a role. These UV photons can be 
important not only in the desorption of ices but also in the cre- 
ation of reactive photo-products such as energetic H atoms and 
radicals which can move through the ice and encounter other 
speci es leading to the fo rmati on of more complex m olecules 
(e.g.. ld'Hendecourt et al.lll982l:lGarrod & Herbstll2006h . 

Water ice is the dominant consitu ent of interstellar ices (e.g., 
iGibb et al.ll2000l: lPontoppidanll2006h with an abundance at least 
three orders of magnitude larger than that of gaseous water 
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in cold clouds (ISnell et alJ bOQOl: iBoonman et"aDl2QQ3l) . Thus, 
evaporation of water ice, even at a low fraction, can signifi- 
cantly aff'ect the gaseous water abundance. Recent models of 
translucent and dense clouds invoke photodesorption of water 
ice in the outer regions to explain the gaseous water emission ob- 
served by the Submillime ter Wave Astronomy Satellite (SWAS) 
('Bergin & Melnickl l2QQ5l) . Photodesorption is also used to inter- 
pret the tentative detections of HDO and oth er gaseous species 
in the surface layers of pr otoplanetary disks dWillacy & Langej 
I2QQQI: IPominik et aDl2005i) . The adopted desorption efficiencies 
in these models, about 0. 1% per incident UV p hoton, are based 
on a single experiment by Westley et al.' ('1995a b) exposing ices 
at 35-100 K to Lyman-a^ radiation. The authors did not detect 
any water photodesorption in the limit of low UV photon flu- 
ence (integrated flux). Therefore it was suggested that water 
photodesorption only occurs from ices that have been subject to 
large doses of UV photons and not directly upon the first expo- 
sure to UV radiation. This is in contrast with recent experiments 
reported by Oberg et al. (submitted to ApJ), where there is a clear 
component of the water photodesorbed from amorphous ices at 
18-100 K that is detected directly upon the first exposure to a 
UV lamp. Apart from that the results remain quite similar to the 
ones by Westley et al. 

Other experi ments on UV irradiation of water ices have also 
been performed. iGhormley & H ochanadel (1971) observed H, 
OH, and H2O2 following UV irr adiation of crystalline ice at 263 
K, while iGerakines et aP (Il996l) found production of OH, HO2, 
and H2O2 in the ice upon exposing amorphous ice at 10 K to UV 
light covering mainl y the first and s econd electronic absorption 
bands of H2O. Watan abe et al.l ("2000) irradiated amorphous D2O 
ice at 12 K with UV photons and observed substantial amounts 
of D2 after irradiation at A = 126 nm , but very lit tle at A = 
172 nm. In the experiments by Yabushita et al. (2006) H atoms 
were found to desorb from the ice after UV irradiation at A = 
157 nm and A = 193 nm. There are also a few reports on two- 
and multi-photon excitat ion of water ice leading to photodesorp- 
tion of H2O molecules (iNishi et al.lll984l: iBergeld & Chakarovl 
2006). In these cases the photon energies are below the thresh- 
old for absorption in the ice, but upon multiple absorptions the 
excitation energies fall between 9 and 10 eV, between the first 
and second absorption band in ice. Clearly, there is a need for 
more quantitative information on the processes induced by UV 
photons in ices, even for the simplest cases such as pure water 
ice. 

We present here the results of the first theoretical study of 
the dissociation of H2O molecules in pure water ice following 
absorption by UV photons. In addition to providing probabilities 
for desorption to be used in astrochemical models, these simula- 
tions provide insight into the mechanisms leading to desorption 
as well as the movements of the energetic photoproducts in the 
ice before they become trapped. In Sect. 2, we present the meth- 
ods used in this study, in Sect. 3 the main results, and in Sect. 4 
a short discussion and astrochemical implications. In Sect. 5 the 
results are summarized and some concluding remarks are given. 



2. Method 

All our calculations have been per formed using classical 
Molecular Dynamics (MD) methods ( Allen & Tildeslevlll987"b 
with analytical potentials. Details of the computa tional proce- 
dure have been described in I Andersson et al.l (l2006h : here only a 
brief outline of the methods will be presented. 



2.1. Amorphous water ice 

To create an amorpho us ice slab, the procedure outlined in 
lAl-Halabi et al.l (l2004ah was used. In brief, a slab of 8 bilayers 
(16 monolayers) of crystalline ice was first created consisting of 
a cell containing 480 H2O molecules. The cell has the dimen- 
sions x: 22.4 A, y: 23.5 A, and z: 29.3 A. Periodic boundary 
conditions are applied in the x- and }7-directions, the z coordi- 
nate being parallel to the surface normal. Thus an infinite ice 
surface is created. The H2O molecules are treated as rigid ro- 
tors and their interacti ons are governed by the TIP4P potential 
( Jorgens en et al.|[T983l) . which describes the interaction as a sum 
of pair interactions (electrostatic and Lennard- Jones potentials). 
The two bottom bilayers are kept fixed to simulate bulk ice and 
the molecules in the other 6 bilayers are allowed to move without 
any dynamical constraints other than that they remain rigid. The 
dynamics are at all times governed by classical Newtonian me- 
chanics. To force the transition to amorphous ice, the surface is 
initially allowed to equilibrate for 5 ps at 10 K but then the tem- 
perature is incr eased to 300 K using a computational equivalent 
of a thermostat (iBerendsen et al .119841) . In this way the top bilay- 
ers form a liquid. The system is left to equilibrate for 100 ps after 
which it is rapidly cooled to 10 K. Then it is once again equili- 
brated for 100 ps. The resulting amorphous ice structure most 
closely resembles the structure of compact amorphous ice ob- 
tained experimentally and is thou ght to be representative of the 
structure of interst ellar v^ater ice (lAl-Halabi e t al 2004a: b). See 
Fig, 4 of Al-Halab i et al.1 (l2004ah and Fig. 2 of Andersson etaP 
(12006h for images. It does not exhibit the m icroporous struc- 
ture that is obtain ed in vapor deposited ice (iMayer & Pletze j 
119861: [Kimmel et a l. 2001a, b). Given the dimensions of the sim- 
ulation cell such a structure is simply not possible to obtain, 
since the pores should have a size on the same order as the cell 
we use. However, we believe that it is a good representation of 
an amorphous ice surface on a local scale, whether that be at 
the "outside" of the surface or inside a void deeper in the ice. 
Implications of this for the obtained results are discussed in Sect. 

m 

In the rest of the paper we will discuss the depth into the 
ice in terms of monolayers. To avoid confusion our definition of 
monolayer is the thickness of ice corresponding to half a crys- 
talline bilayer, i.e., if the ice were crystalline each bilayer would 
consist of two monolayers. For ease of definition the monolay- 
ers have been taken to be divided according to the z values of the 
centers-of-mass of the molecules, e.g., the 30 molecules with the 
largest values of z constitute the top ("first") monolayer. 

2.2. Initial conditions 

Once the ice surface is set up, one H2O molecule is chosen 
to be photodissociated. This molecule is then made completely 
flexible and its intramolecular (internal) interactions are gov- 
erned by an analytic potential energy surface (PES) for the first 
electronically excited state (the A^Bi state) of gas-phase H2O 
based on high-quality ab in itio electronic structure calculations 
(iDob byn & Know leslll997|) . This excited potential is fully repul- 
sive so that absorption into this state leads to dissociation of the 
H2O molecule into H -h OH. The intermolecular interactions of 
the excited state H2O with the surrounding H2O molecules are 
governed by specially devised partial charges for describing the 
electrostatic interactions. In short, a charge of -0.2^ is put on 
the O atom and charges of -hO.l^ on the H atoms. This gives a 
smaller dipole moment than that of the ground state H2O po- 
tential. The eff'ect of this is that a less favorable interaction is 
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obtained with the surrounding H2O molecules, giving higher ex- 
citation energies than with the uncorrected gas-phase potential 
energy surface. This leads to the blue- shift of about 1 eV of the 
ice UV spectra seen in Fig. [TJ which agree well with the first 
UV absorption band in amorphous and crystalline ices. If the 
ground state partial charges were to be used for the excited state 
the excitation spectrum would coincide with the gas-pha se UV 
spectr um. For more details on the potentials see Andersso n et al.l 
(l2QQ6h . Similar procedures using potential energy surfaces for 
higher excited states (B, C, etc.) should in principle give reason- 
able representations of higher-lying absorption bands in water 
ice. 

The initial internal coordinates and momenta of the atoms 
in the selected molecule are sampled by a Monte Carlo proce- 
dure using a semi-classical (Wigner) phase-space distribution 
ISchinke 1993) that has been fitted to the ground- state vibra- 
tional wave function of H2O (Ivan Harrevelt et al.l 120011) . This 
procedure gives initial conditions that are very similar to those 
found in fully quantum mechanical methods and has been shown 
to work well for the description of photodissociation processes 
of gaseous molecules. The transition dipole moment function, 
which governs the strength of th e absorption, is taken from 
the cal culation for gaseous H2O by Ivan Harrevelt & van Hemeri 
(l200Qh . 

Dissociation of molecules in the top six monolayers has been 
considered. For each monolayer all 30 molecules have been dis- 
sociated, one molecule at a time. For each molecule 200 con- 
figurations and momenta were sampled from the Wigner distri- 
bution. This gives 6000 trajectories per monolayer and 36000 
trajectories in total. 

2.3. Calculation of spectra 

The excitation energy is computed by taking the energy diff'er- 
ence between an ice slab with an excited state H2O and one with 
a ground state H2O (with the same coordinates). Each excitation 
is assigned a weight calculated as the square of the coordinate- 
dependent transition moment. By summing the weights of the 
excitation energies binned in 0.05 eV-wide energy intervals, "in- 
tensities" are obtained. Taken together these intensities form a 
UV absorption spectrum for the ice. The monolayers 5 to 6 were 
found to be converged to a "bulk behavior" ( Andersson et al.l 
[2Q05h and could therefore be used to compare the calculated 
spectra to experimental data. The gas-phase spectrum presented 
in Sect. [3] was obtained using the same intramolecular potential 
surfaces as above, but without the surrounding molecules and 
with 1000 sampled configurations. 

2.4. Dynamics of the dissociating molecule 

After putting the molecule in the excited state, the dissociating 
trajectory is integrated with a timestep of 0.02 fs. A maximum 
time of 20 ps has been used before terminating the trajectory. 
Most of the trajectories (99.6%) were terminated before that be- 
cause the system was found in one of the final outcomes (see 
Sect. [3]) with negligible probability of transforming into a diff'er- 
ent state. 

When the excited H2O dissociates, the intermolecular inter- 
actions are smoothly switched into separate interactions between 
the photoproducts and water ice, i.e., H-H2O and OH-H2O po- 
tentials. All details of the potentials and the f unctions used to 
switch between diff'erent potentials are given in I Andersson et aP 
(I2OO6I) . The switching functions connect the partial charges. 



the dispersion interactions and repulsive potentials between the 
H2O-H2O potentials and the OH-H2O and H-H2O potentials af- 
fecting the dissociating molecule. The switches are functions of 
the OH distances (Ron) within this molecule and will give the 
interaction parameters as continuous functions in the range 1 . 1 
- 1.6 A in Ron. The intramolecular potential is switched to the 
ground- state PES, which allows H and OH to recombine to form 
H2O. This switch is done in an analogous way as for the inter- 
molecular interactions above, but here the range of Ron where 
the switch is made is 3.0 - 3.5 A. In this range the excited-state 
and ground-state PES are near-degenerate, so a high transition 
probability between the two states is quite probable. Once Ron 
becomes larger than 3.5 A the system will remain on the ground- 
state PES, even if Ron again becomes smaller than 3.5 A. This 
is what allows for recombination of H and OH. The intermolec- 
ular interactions for the recombined ground state H2O with the 
surrounding H2Q mole cules are taken from the TIP3P potential 
("Jorgensen et al.l ll983h . 

A slightly diff'erent stop criterion has been used in these 
calculations compared to t he results presented previously 
(lAndersson etal.ll2005l l2006h . When an H atom or OH is ac- 
commodated to the ice surface ("trapped") the trajectory is run 
until its translational energy equals T or lower and the bind- 
ing energy to the surface is 0.02 eV (H atom) or 0.1 eV (OH) or 
stronger. In the older version of the code, the stop criterion was 
based on the kinetic energ y of the indiv idual atoms in relation 
to the potential energy ( And ersson et al.ll2 006). The introduction 
of the new termination scheme led to a reduction by about 50% 
of the number of trajectories exceeding 20 ps. 

Although most of the results presented here focus on amor- 
phous ice, calculations have been performed f or crystalline ice 
as well for comparison. Details can be found in lAndersson et al.l 
(i2006D . 

3. Results and Discussion 

3.1. Ice UV absorption spectrum 

As presented in Fig. [H our calculated spectra of the first UV 
absorption bands in amorphous and crystalline ice match very 
well the experimentally obtained spectra, both in general shape 
as well as in the peak and threshold energies. The calculated 
gas-phase spectrum of the first absorption band shown in Fig. 
[T]also matches the experimental peak energy (7.4-7.5 eV) quite 
well. This is naturally to be expected, since the potential surfaces 
used are based on very high-quality ab initio calculations of the 
energy. The success in reproducing the measured spectra leads 
us to believe that the amount of excess energy released into the 
ice is basically correct. 

The ice spectrum is blue- shifted with respect to that of 
gaseous H2O and has significant cross section only in the 7.5- 
9.5 eV range. Thus, the photodesorption probabilities computed 
here are appropriate for photons in the 1300-1500 A range. 
Dissociation of H2O can also occur following absorption into 
higher excited states (e.g., the equivalent of the B state of 
gaseous H2O) but these generally contribute less than 20% of 
the total absorption in a dense cloud. 

3.2. Photoprocesses and desorption probabilities 
3.2.1. Overall probabilities 

Photodissociation of a water ice molecule can have several 
outcomes, with the H and OH photoproducts either becoming 
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Fig. 1. Calculated and experimental (iKobayashJ 19831) spectra of 
the first UV absorption band in crystalline and amorphous ices 
and calculated first absorption band for gas-phase H2O. Note that 
the intensity is given in arbitrary units and that the intensities of 
the ice spectra have been scaled to roughly coincide. 



trapped in the ice, recombining back to an H2O molecule, or 
desorbing from the ice surface. Fig. [2] shows the probabilities 
(as fractions of the number of absorbed photons) of the main 
processes as functions of how deep into the ice the dissociat- 
ing molecule initially is located. Note that these probabilities are 
given per absorbed UV photons and not per incident photon. In 
the first monolayer the dominant outcome is that the hydrogen 
atom desorbs and the OH radical is trapped in or on the ice. The 
probability of this event drops rapidly from about 0.9 in the first 
monolayer to just over 0.1 in the sixth monolayer. This reflects 
the efl'ect of the ice on the motion of the H atom. At the surface 
there are very few molecules to stop the H atom from desorbing, 
but starting from deeper into the ice there are more obstacles on 
the way to the gas phase. 









■ H atom desorbs and OH trapped 
a H + OH trapped 

□ H20 recombined and trapped 

■ H20 desorbs (total) 
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Fig. 2. Fractions of the main outcomes after H2O photodissoci- 
ation for the top six monolayers of amorphous ice. These prob- 
abilities are calculated from all trajectories, irrespective of exci- 
tation energy. 

For the same reasons the probabilities of the other two ma- 
jor outcomes steadily increase as one moves deeper into the ice. 



Table 1. Total probabilities of H atom, OH, and H2O desorption 
(per absorbed UV photon) as functions of monolayer 



ML H atoms OH H2O 



1 


0.92 


0.024 


7.3x10-3 


2 


0.70 


0.015 


7.7x10-3 


3 


0.51 


4.0x10-3 


4.8x10-3 


4 


0.30 


0.00 


3x10-4 


5 


0.21 


0.00 


3x10-4 


6 


0.12 


0.00 


0.00 



These are the events when either both H and OH become sepa- 
rately trapped in the ice or when H and OH recombine to form 
H2O to subsequently become trapped in the ice. Except for the 
top two monolayers the probabilities of these two events are 
roughly equal and in the sixth monolayer the probabilities are up 
to 0.4. The reason for the lower probability of the H2O molecule 
being recombined and trapped in the top monolayers can be un- 
derstood from the open structure of the uppermost layers, which 
more easily allows for the photofragments to escape the region 
of the ice where they were initially formed. 

The probability of photodesorption of H2O is seen to be low 
compared with the above processes, 0.7% in the top layer and 
0.8% in the second layer, and then decreases with distance from 
the surface (see Table [1]). However, this is only considering the 
desorption of intact H2O molecules. If one is interested in the re- 
moval of H2O from the surface without considering what enters 
the gas phase, the dominant mechanism in the top two mono- 
layers is actually desorption of separate H and OH fragments. 
Desorption of OH is in most cases accompanied by the desorp- 
tion of an H atom, but a minor fraction of OH desorption occurs 
with the H atom being trapped in the ice (see Sect. 13.2.2b . In 
summary, the desorption probabilities of OH and H2O are about 
2 orders of magnitude lower than that of H atoms with OH des- 
orption being about twice as probable as H2O desorption if one 
sums over the probabilities from all monolayers (see also Sect. 
13. 

3.2.2. Mechanisms and their probabilities 

Analysis of the trajectories shows that there are three dis- 
tinct mechanisms for H2O removal (see Fig. [3]). Note that 
t hese snapshots a re taken from calculations on crystalline ice 
dAndersson et al.l |2006) for ease of visualization: (a) An H 
atom released from photodissociation of H2O is able to transfer 
enough momentum to one of the other H2O molecules to "kick" 
it off the surface, (b) H and OH recombine to form H2O and 
subsequently desorb, and (c) the H and OH both desorb from the 
surface separately. In Table [2] the absolute probability of H2O 
desorption is given along with the relative probabilities of the 
three different mechanisms for each monolayer. With the addi- 
tional mechanism (c) the probability of H2O removal is 2.7% and 
1.9% per absorbed UV photon in the first and second monolayer, 
respectively. 

In the top two monolayers the direct desorption of H and 
OH fragments is the dominant desorption mechanism, but in the 
third layer the three distinct desorption mechanisms ("a", "b", 
and "c") are roughly equally probable. Following UV absorption 
in monolayers 4 and 5 only the indirect "kick-out" mechanism is 
effective. In Fig.|4]the probabilities of all mechanisms of desorp- 
tion of H atoms, OH radicals, and H2O molecules are presented 
for the top five monolayers. The desorption of OH is possible 
either together with the H atom as shown above or separately 
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amount of matter that has been observed to desorb following 
photoexcitation. For H2O photoexcited below the fifth mono- 
layer there is no evidence of photodesorption of H2O molecules. 



Fig. 4. Fractions of the detailed outcomes after H2O photodisso- 
ciation for the top five monolayers of amorphous ice. The error 
bars correspond to a 95% confidence interval. 
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■ Monolayer 1 
S Monolayer 2 

■ Monolayer 3 
□ Monolayer 4 
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Fig. 3. Snapshots of trajectories of mechanisms of H2O des- 
orption for a crystalline ice model, (a) One of the surrounding 
molecules desorbs. (b) The photofragments H and OH recom- 
bine and desorb as H2O. (c) The photofragments both desorb as 
separate species. The red and white atoms correspond to O and 
H atoms in the surrounding molecules and the blue and yellow 
atoms correspond to the O and H atoms of the photodissociated 
H2O molecule. 

Table 2. Absolute probabilities for removal of H2O from an 
amorphous ice surface (per absorbed UV photon) as functions 
of monolayer 



ML 




(a) 


(b) 


(c) 




desorption 


H2O 


H2O 


H + OH 




probability 


intact 


intact 


fragments 


1 


0.027 


0.10 


0.17 


0.73 


2 


0.019 


0.16 


0.24 


0.60 


3 


7.0x10-3 


0.43 


0.26 


0.31 


4 


3x10-4 


1.00 


0.00 


0.00 


5 


3x10-4 


1.00 


0.00 


0.00 



^Probability for desorption as H2O or H + OH. 



with the H atom remaining trapped. The former mechanism is 
found to have a higher probability. In total, the desorption of 
OH is about a factor of 2 more probable than the desorption of 
H2O. Below the third monolayer the released OH radicals do not 
have sufficient kinetic energy to make it to the top of the surface 
and desorb. The rightmost three categories constitute a further 
division of the category of indirect H2O desorption. Here "H2O 
indirect desorption" only refers to the cases where a molecule is 
kicked out by an H atom, which subsequently remains in the ice. 
The category "H2O desorption induced by recombination" refers 
to the rare occurence where it is the excess energy from the re- 
combination of the H and OH fragments that kicks the molecule 
off' the surface. The case "H + H2O desorb" refers to when the 
H atom kicks the molecule off' the surface and subsequently des- 
orbs itself. This last category dominates the indirect desorption 
in the ffrst two monolayers and also consitutes the maximum 



The division into monolayers in Table |2] refers to where the 
photoexcited H2O is situated. The H2O molecules that actually 
desorb upon being expelled by an H atom or recombining H2O 
molecule all originate in monolayers 1 (84%) and 2 (16%). In 
most of these cases it is not only the transfer of momentum that 
is eff'ective, but also the repulsive interaction from the photoex- 
cited molecule, which most often is in the near vicinity of the 
desorbing molecule. To illustrate this one can consider the low- 
ering in binding energy of the molecule about to be desorbed. 
In monolayer 1 the average binding energy of all molecules is 
0.9 eV (calculated using the TIP4P potential) and of the des- 
orbing molecules prior to excitation it is 0.8 eV. However, the 
photoexcitation lowers the binding energy of these molecules by 
on average 0.3 eV. About 25% of the desorbing molecules do 
not have their binding energy signiffcantly lowered by excita- 
tion, but are kicked out solely by momentum transfer. If these 
molecules are excluded then the binding energy is on average 
lowered by 0.4 eV. In the second monolayer only about 10% of 
the desorbing molecules do not have their binding energy signif- 
icantly lowered. There the average binding energy is 1.1 eV and 
this is lowered by 0.3 eV on average upon photoexcitation. 

There is no sign of molecules being electronically excited 
and then desorbing intact directly, i.e., expelled by the repulsive 
interaction of the excited state molecule with its surroundings. 
The molecules dissociate very quickly (on the order of 10 fs) 
and that is not sufficient time for the molecule to desorb before it 
is dissociated. As discussed above the photofragments can how- 
ever recombine and then subsequently desorb as H2O. 

3.2.3. Effects of product energies 

The H atoms that are released have an average energy of 1.5- 
2.5 eV depending on the excitation energy and to a lesser extent 
in which monolayer they originate (see Fig. [5]). As can be seen 
the average initial H atom translational energy increases with 
increasing excitation energy, but above an excitation energy of 
9 eV it drops to somewhat lower energies. The average vibra- 
tional energy of OH has a minimum value of 0.3 eV around Eq^^c 
= 8 eV but increases strongly to about 2 eV at £^exc = 9.5 eV. 
This implies that at lower excitation energies the vast major- 
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Fig. 5. The average initial translational energies of H atoms 
[^trans,i(H)] and OH radicals [£^trans,i(OH)] and the average OH 
vibrational energy following photodissociation as function of ex- 
citation energy for the top six monolayers in amorphous ice. 



ity of the OH molecules are formed in the vibrational ground 
stat e, since the e xperimental zero-poi nt energy of OH is 0.2 3 
eV (iHuber & He rzberg 1979) (see also I Andersson et al.l (l2006h ). 
When the excitation energy is increased above 9 eV, large frac- 
tions of vibrationally excited OH are produced. The average ini- 
tial translational energy of OH is only weakly dependent on ex- 
citation energy and lies around 0.2 eV with only a slight increase 
with increasing £^exc- 
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Fig. 6. The probabilities of H atom desorption and trapping of 
H and OH as function of initial H atom translational energy av- 
eraged over the top six monolayers in amorphous ice. Note that 
for technical reasons these are probabilities with the possibility 
of H2O recombination excluded (see text). The error bars corre- 
spond to a 95% confidence interval. 

In Fig. [6] the H atom desorption probability is plotted along- 
side the probability of H and OH both becoming trapped as func- 
tions of initial H atom translational energy. The plotted probabil- 
ities are somewhat higher than they should be because the prob- 
ability of recombination of H2O has been excluded in the set of 
outcomes. The reason for this is that during recombination the 
translational energy of the H atom becomes very high. If recom- 



bination occurs immediately after dissociation it is quite diffi- 
cult to distinguish the maximum translational energy the H atom 
normally would have after photodissociation and the maximum 
translational energy it gets during recombination. However, the 
trend is clear that the H atom desorption probability increases 
with increasing initial translational energy, as one would intu- 
itively expect. The reason for the unexpectedly high desorption 
probability at £^trans,i(H) = 0.7 eV is not quite clear and it could 
simply be an effect of insufficient sample size, given that the er- 
ror bars are fairly large. 

Similarly, the desorption probability of OH has been plotted 
in Fig. [7] as function of initial OH translational energy in the top 
three monolayers. Also in this case the desorption probability in- 
creases with increasing initial kinetic energy. The effect is much 
stronger than for the case of H atoms, which reffects the much 
stronger binding energy of OH to its surroundings compared to 
that of the H atom. Interestingly, if the desorption probability 
is weighted with the initial distribution of translational energies, 
the total OH desorption probability (for the top three monolay- 
ers) increases rapidly at 0.2 eV (the average initial translational 
energy) and attains a basically constant value of 0.1% for all en- 
ergies above that. 
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Fig. 7. The probabilities of OH desorption as functions of ini- 
tial OH translational energy averaged over the three top mono- 
layers of amorphous ice. Also shown are the initial distribution 
of OH translational energies and the total desorption probabil- 
ity (desorption probability x initial distrbution). The error bars 
correspond to a 95% confidence interval. 



The dependence of the indirect H2O desorption probability 
on the translational energy of the H atom is found to be weak (see 
Fig. [8]). There is not much evidence of any variation with trans- 
lational energy and the desorption probability is fairly constant 
at 0.1% (averaged over the top six monolayers) over the whole 
energy range. It would be natural to expect that there could be 
a strong dependence on translational energy, but as discussed 
above the desorption is in most cases a combined effect of a re- 
pulsive force from the excited molecule, a lowered binding en- 
ergy, and the momentum transfer from the H atom. Apparently, 
this allows also H atoms with relatively low translational ener- 
gies to kick out H2O molecules. 
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Fig. 8. Probability of indirect H2O desorption ("kick-out") as 
function of initial H atom translation energy (divided into bins of 
0.5 eV intervals) averaged over the top six monolayers of amor- 
phous ice. The error bars correspond to a 95% confidence inter- 
val. 
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Fig. 9. Probabilities of desorption of H atoms, OH radicals, and 
H2O molecules as function of excitation energy averaged over 
the top six monolayers of amorphous ice. The error bars corre- 
spond to a 95% confidence interval. 
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3.2.4. Dependence on photon energy 

Considering the desorption probabilities as functions of the ex- 
citation energy (Fig. [9]) it is interesting to note that H atom des- 
orption becomes less probable with increasing excitation energy 
(0.8 at 7.3 eV and 0.3 at 9.5 eV). Since it was shown that the av- 
erage initial translational energy mainly increases with increas- 
ing excitation energy (Fig. [5]) and that the desorption probability 
increases with increasing translational energy (Fig. [6]) this seems 
like a paradox. The simple explanation of this behavior is that the 
lower excitation energies dominate in the top monolayers while 
the more energetic UV photons are mainly absorbed towards the 
bulk of the ice (see, e.g.. Fig. 5 of Andersson et al. 2005). The 
desorption probability summed over the whole excitation energy 
range decreases rapidly with depth into the ice (Table [B and 
therefore this unexpected behavior is found. The desorption of 
OH seems to increase with increasing excitation energy. This 
is a reflection of the fact that even though the average initial 
translational energy of OH varies only slightly over the excita- 
tion energy range (Fig. [5]), there is a high-energy tail of the OH 
translational energy distribution (see Fig. [7]) that becomes larger 
with higher excitation energies. The desorption probability of 
H2O does not show strong dependence on excitation energy, but 
considering the rather large error bars some energy dependence 
cannot be entirely ruled out. 

3.3. Mobility of photoproducts 

The H atoms produced in the photodissociation event are found 
to be quite mobile in the ice. On average the H atoms that be- 
come trapped move 8 A from their original locations. In extreme 
cases distances over 70 A are recorded. The OH radicals formed 
in the ice move only about 1 A with maximum distances moved 
of 5 A. However, OH radicals formed from photodissociation in 
the top three monolayers are able in some cases to move tens of 
A on top of the surface (up to more than 60 A). The fact that 
some of the photofragments are able to move large distances im- 
plies an increased probability of reactions with other species in 
or on the ice, than if they would remain in the immediate vicin- 



ity of their point of origin. For more details see I Andersson et al] 
(l2QQ6h . 



3.4. Comparison to experiments 




0.5 1 1.5 2 2.5 3 3.5 4 

Kinetic energy [eV] 

Fig. 10. Calculated distributions of H atom initial translational 
energies and desorption energies upon UV absorption at £^exc = 
7.9 eV summed over the top six monolayers in amorphous ice. 

In the experiments by lYabushita et al.l ('2006') on UV irradia- 
tion of poly crystalline and amorphous ices at 100 K absorption at 
A= 157 nm (£exc = 7.9 eV) was found to result in H atom des- 
orption. The translational energy distribution of the desorbing 
atoms was observed to consist of three components at 0.61 eV, 
0.081 eV, and 0.014 eV, respectively. The lowest-energy com- 
ponent seems to consist of atoms that have thermalized prior to 
desorption. Thermal desorption occurs on a time scale that is 
likely much longer than would be feasible to do with molecular 
dynamics simulations. Therefore, it is not likely that we would 
see this third component in our calculations, but the other two 
should be possible to reproduce. In Fig. [TOlthe calculated dis- 
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tributions of initial translational energies and desorption ener- 
gies of the released H atoms following excitation at 7.9 eV are 
shown. It is seen that the initial translational energy has a peak 
at around 1.9 eV and the desorption energy peaks at about the 
same energy. This clearly is much higher than found in the ex- 
periments, so it seems the energy of the desorbing H atoms is 
overestimated in our calculations. There could be three expla- 
nations to this behavior: (i) either the H atoms lose more en- 
ergy prior to desorption or (ii) they are initially formed with less 
translational energy or (iii) a combination of these two effects. 
The possible sources of loss of highly energetic H atoms in the 
ice that cannot be treated by our calculations are: (a) the loss of 
energy by excitation of intramolecular modes in collisions with 
H2O molecules and (b) reactions with H2O m olecules to form, 
e.g., H 2 and OH. For a thorough discussion see 'Andersso n et all 
(l2006h . Since we cannot directly tell how effective mechanism 
(i) is it is hard to speculate how much the desorption energy 
distribution is cooled through energy transfer and reaction. It 
is easier to speculate about mechanism (ii), since it is possible 
to monitor the dependence of the average desorption energy on 
the initial translational energy. If the initial translational energy 
is around 1 eV the average desorption energies lie in the range 
0.4-0.8 eV, which would be in much better accord with the ex- 
perimentally measured desorption energies. If this is the most 
important mechanism for cooling the desorption energy distriu- 
bution, then at £^exc = 7.9 eV the initial H atom translational en- 
ergies are overestimated by roughly 1 eV. As discussed in Sect. 
1 3. 21 this would probably have little importance for the H2O des- 
orption probability (Fig. [8]). However, the desorption probability 
of H atoms could be somewhat lower than predicted in our cal- 
culations, but it is still quite likely to be of the same order of 
magnitude (Fig. [S]). If the average H atom translational energy 
is overestimated it is also likely that the OH translational energy 
is somewhat overestimated. Considering the weak dependence 
on excitation energy (Fig. [5]) the average OH translational en- 
ergy might not be highly overestimated, but the high energy tail 
could be smaller, meaning that less OH radicals desorb than pre- 
dicted here. If the initial kinetic energies of the photofragments 
are overestimated, the only possibility to account for the blue 
shift of the excitation energy is that the intermolecular repulsion 
is underestimated. This could have as an interesting effect that 
the indirect desorption of the surrounding molecules could be 
underestimated, since they would experience an even larger re- 
pulsive force from the excited molecule than predicted by our 
calculations. This remains to be investigated. 

The high probabilities of H atom desorpti on is also in ac- 
cordance with the finding of Gerakines et all ((1996^ that their 
UV-irradiated water ice was most likely depleted of H atoms, 
since a large amount of oxygen rich products was found. Our 
results on mobility o f the released photofragments (Sect. 3.3, 
lAndersson et al.' 2006) have been supported by recent experi- 
ments by Files et al. ( 2007), who observed an average separation 
of H and OH of 7+2 A immediately following photodissociation 
of H2O upon UV irradiation of liquid water. This is in excel- 
lent agreement with our calculated value of 8 A for the average 
distance from the site of photodissociation of the H atoms in 
amorphous ice. Even though the temperatures are quite different 
in the two cases, it is expected that liquid water and compact 
amorphous ice are quite similar on the short time scales during 
which photodissociation takes place. 

A di rect comparison with the water ice photodesorption re- 
sults by IWestley et all (Il995a|[9) is difficult to make since they 
used Lyman-Qf radiation, which leads to excitation to a higher 



absorption band than considered here (iKobay ashil 1 1 98 3l) . Their 
findings of basically no photodesorption of intact H2O at low 
temperatures in the limit of single-photon absorption can there- 
fore be neither refuted nor confirmed by our results. However, 
the detection of desorbing H2 and O2, and pos sibly OH and 
H2O2 actually agree with the present results (IWestley et al.l 
1995a) (see below). In the experiments of Oberg et al. (submit- 
ted to ApJ) the photodesorbing material is detected in the form 
of OH, H2O, H2, and O2. This is the first time a positive detec- 
tion of OH photodesorption has been reported. The detection of 
OH agrees nicely with our simulations (see also Sect.©. 

It is important to bear in mind that even though only H, OH, 
and H2O would desorb upon absorption of a single UV photon, 
in the experimental setup a large amount of UV photons may 
be absorbed in the ice during a relatively short time interval. 
This makes it possible to produce significant amounts of H and 
OH in the ice, and possibly also O atoms from for instance the 
photodissociation of the formed OH radicals. Given that ther- 
mal diffusion is effective in moving these reactive species into 
close contact, all the aforementioned desorbing species can be 
accounted for through recombination followed by desorption. 

4. Discussion and astrophysical implications 

Based on the results presented in this paper some important con- 
clusions can be drawn on the possible outcomes of UV irradi- 
ation of water ice in interstellar environments. First, it is im- 
portant to realize that by far the most likely species to desorb 
are H atoms, followed by OH radicals. In addition, when there 
is removal of H2O from the surface it seems likely that most 
of it comes off in the form of separate H and OH fragments. 
Therefore, there is not a one-to-one correspondence between 
H2O molecules removed from a surface and H2O appearing in 
the gas phase. 

So far, this paper has been concerned with probabilities of 
photoinduced processes following absorption of one UV photon 
in a specific layer in the ice. However, not all incident UV pho- 
tons are absorbed by molecules in the top six monolayers. To 
estimate desorption probabilities per incident photon rather than 
per absorbed photon one needs information about the absorp- 
tion cross section. In our semiclassical simulations we are not 
able to calcu l ate the absolute absorption cross section. However, 
Mason et al.l (l200tf ) have measured the absorption cross section 
of water ice at 25 K and found the peak absorption cross section 
in the first absorption band to be about 6x10"^^ cm^ at an exci- 
tation energy of 8.61 eV. This leads to an absorption probability 
of 0.007 photons ML"^ (see Appendix lAl for an outline of this 
calculation). 

From the above estimate of the absorption probability it is 
possible to calculate photodesorption probabilities per incident 
UV photon. This is done by weighting the desorption prob- 
abilities per absorbed photon for each monolayer by the ab- 
sorption probability for the specific monolayer with the absorp- 
tion probabilities in any upper monolayers subtracted from the 
incoming photon flux. Using the information in Table [2] one 
arrives at a probability of removal of H2O from the ice of 
3.7 X lO""^ photon"^ About 60% of the removed H2O comes off 
in the form of H -h OH, 20% desorbs as recombined H2O, and 
the remaining 20% consists of H2O "kicked" out from the sur- 
face. The total photodesorption yield of intact H2O molecules 
is 1.4 X 10""^ photon" ^ For OH desorption the probability is 
3x10"^ photon" \ which includes both desorption with and with- 
out H atoms. The H atom photodesorption probability is rela- 
tively high, 0.02 photon" ^ The above ratios of photodesorption 
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of OH and H2O are in excellent agreement with the recent ex- 
periments by Oberg et al. (submitted to ApJ), which inferred 
that roughly equal amounts of OH and H2O photodesorb at low 
surface temperature (18 K). They obtain a total photodesorption 
yield of about 1.3 x 10"^ in the low-temperature limit, which 
is about 3 times higher than what is found from our simula- 
tions. Given the experimental uncertainties and the approxima- 
tions made in the simulations this can be considered as good 
agreement. Commonly used estimates of H2O photodesorption 
probabilities in the range 1 x 10"'^-3.5 x 10"^ have been used 
to model different environments in agreemen t with observations 
(lBerginetal ."1995': "Willacv & Langei^ '2000; S nell et al.1 [20051: 
iDominik et al. 2005; Bergin & Melnick 2005). Our results in- 
dicate that these estimates are reasonable. However, the finding 
that less than half of the desorbed material leaves the grain in 
the form of intact H2O molecules is something that should be 
included in the models. 

The results presented here are only for photoinduced pro- 
cesses that are followed until the photoproducts desorb or are 
thermalized within the ice. This is all happening on a picosec- 
ond time scale. For longer time scales there is the possibility of 
thermal desorption of especially the H atoms. These are rela- 
tively weakly bound to the ice surface and it is quite probable 
that some fraction of the released H atoms desorb after thermal- 
ization in the ice. This contribution to the desorption probability 
is therefore not included in the above estimate. 

It is interesting to discuss the possible eff'ects of the over- 
all morphology of the ice surface. It is clear that the vast 
majority of water ice surface s in the interstellar medium are 
amorphous (Hagen et al.lll98ll) . However, it is debated whether 
the ice is mos tly porous, as found i n vapo r deposited ice 
(iMayer & Pletzer 1986; Kimmel et al. 2001allbh . or more com- 
pact (iFraser et al.h2004i: iPalumba fc2006). A complicating factor 
is that most likely the ice is formed through chemical reactions 
on gra ins rather than accretion o f H2O molecules from the gas 
phase do 'Neill&Williamslll999l) . Therefore, the exact ice mor- 
phology that results from such a chemical b uild-up of the ice is 
not yet clear (see however'Cuppen & HerbstI (12007 )). If a porous 
ice surface is subjected to UV irradiation one could have release 
of H, OH, or H2O into a void in the ice rather than directly into 
the gas phase. This is inferred to happen in the experiments by 
lYabushita et al.l (l2006l) where a large component of the desorb- 
ing H atoms released after UV irradiaton of amorphous ice are 
thermalized, likely due to H atoms being accommodated within 
a void in the ice and then desorbing thermally. Some OH and 
H2O could photodesorb in a similar way. However, if the path to 
reach the gas phase from inside a pore is restricted, these species 
have a high probability of being trapped inside a pore because of 
their strong attractive interactions with the ice. In practice, the 
OH and H2O released in this way might not show up in desorp- 
tion. Therefore, if a porous ice is considered it is necessary to 
distinguish between the surface area that is directly exposed to 
the gas phase and that which is within a void with restricted ac- 
cess to the outside. In conclusion, the H atom desorption from 
a porous ice is likely to be diff'erent from that of compact non- 
porous ice, but the desorption yields of OH and H2O are not 
necessarily very diff'erent in the two types of ice. 

The cosm ic ray induced UV flux inside dense clouds is about 
10 cm" -h-^ (IShen et al.ll200i ). For an ice-coated grain with a 
typical size of 0.1 /im, this would give an arrival rate of about 
1 UV photon per day. The case of a single UV absorption event 
as described in our simulations therefore gives a realistic picture 
of photodesorption in dense clouds. In the laboratory, the UV 
flux is many orders of magnitude higher and multiple absorp- 



tions of UV photons within the ice surface in a short time inter- 
val may drive secondary reactions of photodissociation products 
from difl'erent H2O molecules. Experiments with difl'erent UV 
flux levels down to low levels will be needed to provide quanti- 
tative data on photodesorption yields relevant for astrophysical 
applications. 

As has also been noted in our previous work another impor- 
tant aspect is the release of reactive species into the ice. That 
would have implications for reactivity in the ice with, e.g., CO 
that could react with energetic H or OH to form HCO or CO2. 
This could be one clue to unraveling the mystery of CH3OH 
and CO2 formation in the interstellar medium. Indeed, formation 
of CO2 is readily observed when a mixed H20:C0 ice is pho- 
tolysed (jd 'HendecourtetalJ 119861: IWatanabe & Kouchil l2002t 
IWatanabe et al.ii20Q7E 



5. Conclusions 

We have shown that it is possible to have H2O photodesorption 
upon UV absorption to the first absorption band in the top five 
monolayers of an amorphous ice surface. The main mechanisms 
for this photodesorption are either photodissociation followed 
by recombination of H and OH and subsequent desorption of 
the recombined H2O molecule or a "kick-out" of another H2O 
molecule in the ice by the energetic H atom released from pho- 
todissociation or, less likely, by the energy released from a re- 
combined H2O molecule. In most cases, however, removal of an 
H2O molecule from the ice is in the form of separate H and OH 
fragments. An estimate of the photodesorption yield per inci- 
dent UV photon from our calculations agrees well with the H2O 
photodesorption yields that are commonly used in modeling as- 
trophysical environments. 

UV absorption leads in most cases to desorption of H atoms 
or the trapping of H and OH in the ice either as separate frag- 
ments or as recombined H2O. The desorption of H atoms is about 
2 or 3 orders of magnitude more probable than desorption of OH 
and H2O. The OH desorption probability is about twice the H2O 
desorption probability. The high mobility of H atoms inside the 
ice and OH radicals on the ice surface will facilitate formation 
of other molecules such as CO2. 
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Appendix A: Estimate of absorption probability 

To calculate an absorption probability per monolayer in an in- 
terstellar ice surface one needs an estimate of the effective area 
taken up by one molecule. Since the angle of incidence of the 
photon is arbitrary, all possible incidence angles have to be taken 
into account, not only normal incidence. The surface area of 
our simulation cell is 22.4 A x 23.5 A = 526 A^ = 
5.26 X 10"^^ cm^. If one assumes a flat surface, the average 
efl'ective surface area seen by a photon from any incidence angle 
between 0° (normal incidence) and 90° (parallel to the surface) 
is the actual surface area divided by two. This is arrived upon by 
the following expression: 



A cose sin Ode 



.;r/2 



sinede 



A 
2' 



(A.l) 



with the efl'ective area given by Aq^ = A cos e where A is 
the actual surface area and e is the angle of incidence. The av- 
erage efl'ective surface area of the simulation cell is 2.63x10"^"^ 
cm^. Each monolayer consists of 30 H2O molecules. A single 
molecule will therefore have an average efl'ective area (A^^) = 
8.77 X 10"^^ cm^. iMason et al] (l2006l) measured the peak ab- 
sorption cross section, cr, in the first absorption band to be 
about 6x10"^^ cm^ around an excitation energy of 8.61 eV. With 
this value the absorption probability per monolayer becomes 
= (t/(A^^) = 7 X 10"^ (this number may vary somewhat 
with excitation energy). If one considers an infinitely deep ice 
surface all photons will be absorbed within the ice. For an ice 
surface in the interstellar medium consisting of a finite number 
of layers, the value of may vary depending on the shape of 
the grain, since the surface is not necessarily flat. However, it is 
likely that any photon that passes through an ice mantle will be 
absorbed by the silicate grain, leading to almost total absorption 
of the incident photons within the grain. 



